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O ■ Abstract: We review recent results from the HI and ZEUS experiments at HERA on 

. charm and beauty production in ep collisions at 300 - 318 GeV centre-of-mass energy. 



1. Introduction 



> 
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. Deep inelastic scattering (DIS) at HERA offers unique opportunities to test and refine our 

understanding of heavy quark production in terms of perturbative QCD. The dominant 
mechanism here is boson-gluon fusion (BGF): a photon coupling to the scattered positron 
^ I interacts with a gluon from the proton to form a quark- antiquark pair. A quantitative 

V ' description of this process requires the gluon momentum distribution in the proton, a 

partonic matrix element and a fragmentation function. The gluon density is known to an 
• accuracy of a few percent from the analyses of scaling violations of the proton structure 

function F2 measured at HERA The masses of the charm and, even more so, of the 
^ , beauty quark ensure that a hard scale is present that renders QCD perturbation theory to 

^ I be applicable to the calculation of the hard subprocess. Fragmentation functions, which 

account for the long-range effects binding the heavy quarks in observable hadrons, are 
extracted from e~^e~ annihilation data, where the kinematics of the hard process is well 
determined; results with high precision appeared recently [^. Compared with the clean 
e^e~ case, the complication in ep collisions lies in the strongly interacting initial state. 
However, relative to other production environments like hadron-hadron collisions or two- 
photon interactions, uncertainties related to hadronic structure are reduced to a minimum. 

QCD calculations have been performed up to fixed order in the so-called massive 
scheme, where only gluons and light quarks are active partons in the initial state. They 
are available in the form of Monte Carlo integration programs [^], which, by using Peterson 
fragmentation functions @, provide differential hadronic cross sections. Due to the higher 
quark mass, the QCD predictions are expected to be more reliable for beauty than for 
charm. However, we note that the NLO corrections to the predicted DIS cross section are 
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around 40% of the LO result in both cases. At very high momentum transfers, a treatment 
in terms of heavy quark densities in the proton may be more adequate; but differences 
between these schemes are not yet significant in the range covered by HERA so far j^. 



in the NLO DGLAP scheme 
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Figure 1: Charm contribution to the proton 
structure function, compared with NLO QCD. 



2. Charm 

Most of the HERA results on charm make 
use of the "golden" decay channel D*^ — > 
D°TT+ followed by D° K--K+; ZEUS also 
uses semileptonic decays. The contribution 
of charmed final states to DIS is quantified 
as the ratio F2/F2, where is defined in 
an analogous way to the proton structure 
function F2 by 

X, y, and are the standard DIS scaling 
variables. Measurements of F2 by both ex- 
periments ||6|, |7| and by using different chan- 
nels yield consistent results. The charm con- 
tribution F2/F2 is found to be about 20 to 
30 % in most of the kinematic region at HERA. 
It is large where gluon-induced reactions dom- 
inate, and decreases only as becomes smaller than ~10 GeV'^, or at higher x values 
^ 0.01, where the quark content in the proton takes over. Figure || displays the measured 
values of F2. The dependence, for fixed values of x, is steeper than for the inclusive 
structure function. The NLO QCD calculation |p, with a gluon distribution extracted 
from HI F2 data, agrees well with the data, which demonstrates the overall consistency 
of the boson-gluon-fusion picture. At low x, the data tend to be somewhat higher and to 
vary stronger with than the prediction. 

o ^ 1.5 < < 2.5 2.5 < p-r < 4 4<p-r<10GeV 

The available statistics makes more de- 
tailed investigations possible. It was observed 
earlier that the NLO QCD calculations with 
Peterson fragmentation (HVQDIS) do not re- 
produce the rapidity (r/) distribution of the 
produced charm meson well in the forward 
region (the outgoing proton direction) ^ . 
The double differential D* cross section 
displayed in Figure ^ reveals that in the HI 
data this is predominantly a feature of the 

low pt{D*) region. The measurement is also compared with the CASCADE Monte Carlo 
program 1^], based on the CCFM evolution equation ||l^, which resums higher order con- 
tributions at low X. Using an unintegrated gluon distribution extracted from inclusive HI 



□ CASCADE 

□ HVQDIS 



H1 



Figure 2: Double differential D*^ cross sec 
tion, compared with QCD predictions. 
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data, it reproduces the data in the low pT region well, but overshoots at higher pT- Such 
shape differences imply that the extrapolated result for F| is model-dependent, but also the 
F2 prediction depends on the evolution scheme. HI has performed a consistent extraction 
and comparison in the CCFM scheme and found somewhat better agreement in the low x 
region than in the standard Altarelli-Parisi scheme. ZEUS 

A possibility to include higher order processes 5 1.2 
in the modeling of heavy quark production is to | 1 
use the concept of photon structure also at non- ^ o.a 
zero virtuality. One can classify events as "direct" "-^ 

o 

or "resolved" according to the measured value of ^ 0.4 

x'^^^ = {E - Pz)2 jets/{E - Pz)all hadrons in di- i "'^ 

jet events. Keeping the LO photoproduction Ian- ° 
guage, this corresponds to the momentum fraction 
of the incoming parton in the photon, but more 
generally, x^^^ is sensitive to any kind of non- 
collinear radiation in the event. The ratio of resolved vs. direct cross sections has been 
determined in this approach by ZEUS [0] and is displayed as a function of virtuality in 
Figure ^. In contrast to the situation for light quarks, the ratio in the DIS regime is very 
similar to that at ^ 0; as expected, using e.g. the virtual photon structure function 
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Figure 3: Ratio of low to high 
cross section, compared with predictions. 



set SaSlD |12| implemented in the HERWIG program |]J. The CASCADE model, with 
gluon emissions ordered in angle rather than in kx, effectively incorporates the perturba- 
tive, anomalous part of the photon structure and reproduces the data well at all Q^, but 



not the AROMA Monte Carlo program |14], which does not include such contributions. 

Charm production has also been measured sep- 
arately in e~p and e+p DIS by ZEUS [^]. The data 
are shown in Figure ^ as a function of . The e~p 
and e'^p results are only barely consistent with each 
other; for > 20 GeV^, the discrepancy amounts 
to 3 standard deviations. However, both measure- 
ments are compatible with the theoretical expecta- 
tion, in which no mechanism exists to generate an 
asymmetry with respect to the lepton beam charge 
at such low four-momentum transfers. 

In summary, the BGF concept at NLO works 
well for charm in DIS, up to high Q^. The HERA 
data reach the precision to identify regions (e.g. at 
low x), where refinements are becoming necessary. 
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3. Beauty 



Figure 4: D*^ cross section in e p 
and e+p DIS compared with NLO QCD. 



Beauty production at HERA is suppressed by two orders of magnitude with respect to 
charm. All HERA measurements of b production so far rely on inclusive semi-leptonic 
decays, using identified muons or electrons in dijet events. Two observables have been 
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used to discriminate the b signal from background sources. The high mass of the b quark 
gives rise to large transverse momenta pJ^'' of the lepton relative to the direction of the 
associated jet. Using this method, both collaborations have published photoproduction 
cross section measurements p^ , 17 1, which are higher than NLO QCD expectations. 

More recently, with the precision offered by the 10 ^ 
HI vertex detector it has become possible to 
observe tracks from secondary b vertices and to ex- 
ploit the long lifetime as a 6 tag, using e.g. the im- 
pact parameter 6. This improves the photoproduc- 



E 

zi. 
O 

°10 



H1 




Data (prel.) 
Fit 

beauty 
charm + fake 



-0.1 



-0.05 0.05 0.1 
impact parameter 5 [ cm ] 



0.15 



Figure 5: Muon impact parameter dis- 
tribution, with decomposition from the 
likelihood fit. 



tion result [19| and provides a first measurement in 
DIS , where resolved contributions involving the 
non-perturbative hadronic structure of the photon 
are expected to be suppressed The DIS case is 
therefore complementary and theoretically simpler. 
The sensitivity to determine the beauty component 
is maximized by combining both variables in a like- 
lihood fit to the two-dimensional distribution in 5 
and p^K The consistency of the two observables has 
been established with the larger statistics available 
in the photoproduction regime The 5 distribution 

for muons in dijet DIS events selected from a dataset corresponding to 10.5 pb^^, is shown 
in Figure |5| together with the decomposition from the two-dimensional fit, which yields a 
bb fraction of fb = (43 ± 8) %. 

A DIS cross section of cr^** „ = 39 ± 

8 lb 10 pb is extracted in the kinematic range 
given by 2 < Q2 < ^qO GeV^, 0.05 < y < 
0.7, prifi) > 2 GeV and 35° < 0{fj,) < 130°. 
This can be directly compared to NLO QCD cal- 
culations implemented in the HVQDIS pro- 
gram, after folding the predicted b hadron dis- 
tributions with a decay lepton spectrum. The 
result, 11 lb 2 pb, is much lower than the HI 
measurements. The data have also been com- 
pared with the CASCADE Monte Carlo simula- 
tion; the result of 15 pb also falls considerably 
below the measurements. 

We summarize the HERA b results (IT 
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19, 20 1 as a function of Q in Figure g, where 



Figure 6: Ratio of measured b produc- 
tion cross sections at HERA to theoretical 
expectation, as a function of Q^. 



the ratio of the measured cross sections to the- 
oretical expectations based on the NLO QCD 
calculations ||3|, ^ is displayed. It is consistent 

with being independent of Q^. The discrepancy between data and theory is similar to the 
situation observed in pp and, more recently, 77 interactions p^]. The first measurement 
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in DIS indicates that in ep collisions this is not a feature of hadron-hadron like scattering 
alone. 
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